Nonlinear Aharonov-Bohm scattering by optical vortices 
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We study linear and nonlinear wave scattering by an optical vortex in a self-defocusing nonlinear 
Kerr medium. In the linear case, we find a splitting of a plane-wave front at the vortex proportional 
to its circulation, similar to what occurs in the scattered wave of electrons for the Aharonov-Bohm 
effect. For larger wave amplitudes, we study analytically and numerically the scattering of a dark- 
soliton stripe (a nonlinear analog of a small-amplitude wavepacket) by a vortex and observe a 
significant asymmetry of the scattered wave. Subsequently, a wavefront splitting of the scattered 
wave develops into transverse rnodulational instability, "unzipping" the stripe into trains of vortices 
with opposite charges. 

PACS numbers: 42.65.Tg, 03.75.Fi 
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It is known that the magnetic vector potential influ- 
ences the dynamics of a charged quantum particle even 
if the magnetic field vanishes, e.g. when the field is con- 
fined to a cylinder into which electrons can not penetrate. 
This phenomenon is known these days as the Aharonov- 
Bohm effect pi , and it was shown to have a more general 
geometrical framework being directly linked to the wave 
front dislocations and geometrical phases g] . It was also 
shown to have a classical analog in the scattering of a 
linear waves by a vortex |J . Recently, the study of the 
interaction of a water wave with a vortex |J] revealed both 
similarities and differences with the Aharonov-Bohm ef- 
fect, allowing to observe directly the macroscopic aspects 
of the geometrical phases. 

In this Letter we study an optical analog of the 
Aharonov-Bohm effect - the wave scattering by an op- 
tical vortex in a bulk dielectric nonlinear medium. We 
describe, for the first time to our knowledge, nonlinear 
effects in the Aharonov-Bohm scattering including both 
the crossover regime and strongly nonlinear case. Since 
our basic equations describe also the dynamics of a su- 
perfluid flow in the Bose-Einstein condensates M , we ex- 
pect that our results will be useful for the experimental 
study of the similar effects in the nonlinear matter-wave 
physics. 

Model. We consider the propagation of a monochro- 
matic scalar electric field in a bulk nonlinear Kerr-like 
medium with an intensity-dependent refractive index, 
n = hq + 77.2/, where n-i is the Kerr-effect coefficient 
characterizing a defocusing nonlincarity (712 < 0). In 
the paraxial approximation, the Maxwell equations can 
be reduced to the nonlinear Schrodinger (NLS) equation 
for the renormalized slowly varying amplitude ip, 
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coordinates (x,y) are measured in the units of fc |n2|/o 
and (no|n 2 |fc 2 /o) 1 / 2 , respectively. Here fco is the free- 
space wave number, no is the linear refractive index, Vj_ 
is the transverse gradient operator, and / = \E\ 2 is the 
field intensity (see details in Ref. ||). 

Equation (pi) has a stationary solution for a vortex B , 
ipo(r) = Ro(r)e llp , where the real amplitude Ro(r) satis- 
fies the conditions R Q (0) = 0, and R (r) = 1 + 0(l/r 2 ), 
as r — ■> 00. To study the wave scattering by the vor- 
tex, we introduce a new function $(r, z), substituting 
ijj(r,z) = ipQ(r)<&{r,z) into Eq. (|l|) and taking into ac- 
count the lowest-order terms in r~ 1 . The corresponding 
equation for <I>(r, z) has the form, 
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where the propagation coordinate (z) and the transverse 



and it describes the advection of the field $(r, z) by the 
vortex phase gradient field e v /r. 

Linear limit: the Aharonov-Bohm scattering. First, 
we analyze the scattering of a small-amplitude wave on a 
vortex. This problem was studied earlier in the context of 
the superfluid hydrodynamics fij]. Herewith, we develop a 
different approach that allows a simple further extension 
to the nonlinear regime. 

Let us consider a plane wavefront with the coordinate 
x = X(y) that moves in the x-direction with the velocity 
v [see Fig. 1(a)], and is advected simultaneously by the 
vortex velocity field u'°'(jf,z) = —(y/r 2 ) = —y/(y 2 + 
v 2 z 2 ). If the wavefront motion starts at z = z- m < 
and the point z — corresponds to the moment when 
the unperturbed line crosses the vortex [see Fig. 1(c)], 
the shape of the line in the y-direction, h^°'(y, z), can be 
found as (for \y\ ^> 1): 

/ l (o) ( ,,,) = -^Mftan-^-tan-^). (2) 
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Thus, the total deformation of the wave front can be 
calculated as follows (for t>|z|, v|zi n | ^> \y\ ^> 1): 



where 



fc (o) (y , z) = _^M ff . 



(3) 



Equation (g) coincides with the result given by Eq. (4.73) 
of Ref. @ provided v has a meaning of the sound speed. 
It describes a splitting of the wavefront by the amount 
Ax = —(2tt/v), predicted by the Aharonov-Bohm effect 
for the linear waves scattered by a vortex. 
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FIG. 1. Sketch of the wavefront deformation and splitting 
at different distances from the vortex. 

Nonlinear regime: Deformation of a dark-soliton 
stripe. In the nonlinear regime of the vortex-wave scat- 
tering, the steady-state propagation of the wave packet 
in the presence of nonlinearity is possible in the form of a 
solitary wave when nonlinear effects are balanced by the 
wave dispersion |8| . In the problem under consideration, 
such a solitary wave corresponds to a localized intensity 
dip on a constant plane- wave background, a spatial dark- 
soliton stripe ra. The corresponding nonlinear scattering 
problem is rather complicated, and we restrict our study 
by the region of large distances from the vortex core, 
r = 0(e _1 ), where e 4C 1 is a formal small parameter of 
the asymptotic analysis we develop below. 

We assume that, under the action of the vortex field, a 
dark-soliton stripe moves with the velocity v = O(l) and 
the coordinate of its center gets shifted, in respect to its 
unperturbed motion x = vz, by the amount h(ey, ez) = 
0(1). Introducing new variables Y = ey, Z = ez, £ = 
x— vz— h(Y, Z), we look for the solution in the asymptotic 
form, <I> = X^^Lo s n ^ n \ In the zero order in e, we find 



19 2 $ (0) <9$ (0) 

— 771 

2 &e os, 



(1_|$(°)| 2 )$(°)=0, 



that has a well-known solution for a dark soliton 



$ (0) = [VI - v 2 tanh( v / l - v 2 £) + iv]e lt 
In the first order in e, we obtain the equation 
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Y 2 + v 2 Z 2 ' 
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A solution of Eq. (|5|) can be found as $W = 
V{Y, Z)(d$M/dv). Result (|), being expressed as h z = 
V — Y(Y 2 + v 2 Z 2 )~ 1 , has a simple physical meaning: a 
local change of the stripe coordinate, hz, is a sum of 
its local velocity, V, and a correction produced by the 
vortex. 

In the second order in e, we obtain an equation for $^ 2 ' 
similar to Eq. (pi) with different r.h.s. terms. The equa- 
tion governing the evolution of the dark-soliton-stripe pa- 
rameters can be obtained from the orthogonality condi- 
tion of the r.h.s. of the equation for t&W to the solution 
of the homogeneous problem, $W — <9<j>(°) jQ£ m Taking 
into account the parity properties and assuming the nor- 
malization lm(d<&( ' / dv) — 1, we derive the equation 
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where X 2 (v) — |(1— v 2 ). System of Eqs. (||), M) describes 
the evolution of the stripe parameters in the presence of 
a vortex field, and it can be rewritten in the original 
variables (y, z), if we define u = eV, and hz = u + u < -°>, 
as follows 



d 2 u , 9 <9 2 u 
dz 2 dy 2 
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Nonlinear wavefront splitting. The Cauchy problem for 
the elliptic equation (||) is ill-posed, hence its solutions 
are meaningful only if the generation of short- wavelength 
harmonics is not essential. To solve Eq. (g), first we apply 
the Fourier transformation and obtain the equation 



d 2 u k 
dz 2 
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where Uk(z) is the Fourier transform of the field u(y, z). 
Except the special case A = v, the solution of Eq. (M can 
be found in a general form 



Uk 
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(10) 



where a k — iir\ 2 (v 2 — A 2 ) _1 sign(fc), and the values of the 
amplitudes A k and B k are determined from the conti- 
nuity conditions, 



At-A- k =-{Bt-B-) = 



t\v 



{v 2 -\ 2 y 



First, we consider the region z < 0. If the stripe in- 
stability is weak, i.e. A < v, the third term in Eq. jlC| ) 
dominates, and we find the result: u(y,z) w [A 2 /(w 2 — 
\ 2 )]u^(y, z). Thus, the deformation of a dark-soliton 



stripe due to its interaction with the vortex field can be 
presented in the form, 



h(y,z) 



{v 2 - A 2 ) 



h<®( V ,z), z<0, 
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where h^ (y, z) is defined by Eq. (g) of the linear the- 
ory. Therefore, the result ([n]) describes the so-called 
nonlinear Aharonov-Bohm scattering on the vortex; it is 
stronger than its linear analog by the factor v 2 /(v 2 — A 2 ), 
provided A < v as v > 1/2. 

The general solution (FLO) describes a growth of insta- 
bility of a dark-soliton stripe. If this instability is rather 
strong (i.e. A > v), the term with A^ becomes impor- 
tant, and the amplitude coefficients Aj~ and B^ should 
be determined from the initial conditions at z = z m . If 
h(y, z- m ) = h z (y, Z- m ) = 0, in the presence of a vortex the 
function Uk grows linearly, for small z-\-z- ln , and exponen- 
tially, for large z + z m . This type of the stripe evolution 
was earlier discussed in Ref. || for the special case v = 0. 

More importantly, we can find an exact non-singular 
analytical solution of Eq. (|10|) for any v and A. Assuming 
z = Zj n < and selecting the initial conditions in the 
form u(y, z in ) = u z (y, z in ) = 0, h(y, z in ) = 0, h z (y, z in ) = 
u^(y,Zi n ), we solve Eq. (Q) and find the dark-soliton 
deformation as following 
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Exact analytical solution (|12|) is more general than 
the result (^lj), and it transforms into Eq. ( |ll| ) for 
\(v — X)z in /y\ ~^> 1. Solution ([12]) remains valid even for 
v — ► 0, i.e. for deformations of a black soliton. 

In the region z > 0, the exponentially growing term 
~ AJ. of Eq. (nfl) dominates, and the corresponding di- 
vergence at large k is a direct manifestation of the ill- 
posed problem rtq). It can be regularized by the higher- 
order term ~ [3h zyy , which appears in the equation 



h zz +X 2 h 



,(°) 



yy 



due to the effect of radiation taken into 



account in the asymptotic expansion (see details of such 
a technique in Ref. p0|). This term makes the linear 
problem (H) well-posed, but it appears with additional 
nonlinear corrections that all describe the dynamics of 
the transverse instability of the dark-soliton stripe in the 
presence of a vortex field. The asymptotic analysis be- 
comes very complicated and, in order to analyze the long- 
term dynamics of the nonlinear wave scattering, bellow 
we carry out numerical simulations. 

Numerical simulations have been conducted by means 
of a standard beam propagation method, solving Eq. (Q) 
over a numerical grid of 1024 x 1024 points. As initial 



conditions in the nonlinear problem, we use a superposi- 
tion of a vortex and a dark-soliton stripe: tp(x,y,0) = 
ipo(r)®(x,0), where ipo(r) = Ro(r)e i(p and Ro(r) de- 
scribes the stationary solution for an optical vortex soli- 
ton found numerically by a shooting technique. The func- 
tion $(x,0) describes a dark-soliton stripe with an ini- 
tial offset from the position of the vortex (xo), and it is 
given by Eq. (Q) . To ensure unperturbed propagation of 
an optical beam over long distances, the initial function 
V>(x, y, 0) was superimposed on a super-Gaussian-beam 
background with HWHM more than 40 times larger than 
the characteristic scale of the localized structures and in- 
teraction domain. 
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FIG. 2. (a) Snapshots of the interaction of an optical vortex 
with a dark-soliton stripe (v 2 — 0.9). Initially, the stripe is 
shifted at Xo = 7.59 from the vortex and, if propagating alone, 
it reaches the origin at z — 8. For better visualization, 16% 
of the computational domain are shown, and the contrast of 
the images was adjusted, (b) Asymmetric stripe deformation 
at different propagation distances, (c) Vortex horizontal (A x ) 
and vertical (A y ) shift from its initial position. 

In our simulations, we vary the contrast of the dark- 
soliton stripe (i.e. its transverse velocity), and moni- 
tor the stripe deformation for different propagation dis- 
tances. First of all, we investigate the scattering of a 
small-amplitude dark-soliton stripe by an optical vortex, 
which should correspond to the scattering of a linear 
wavepacket associated with the Aharonov-Bohm effect. 
In Fig. 2(a), we show the snapshots of the scattering 
process for the dark-soliton stripe with the transverse 
velocity v 2 — 0.9, at different propagation distances z 
(marked on the figure). Below in Fig. 2(b,c), we show 
the details of the stripe deformation (first symmetric but 
then asymmetric) during the scattering, and the vortex 
shift A induced by the stripe scattering lljfl . 

To investigate the strong nonlinear regime of the 
vortex-stripe scattering we consider different values of the 
stripe amplitude by changing its velocity and, therefore, 
contrast. The results of these calculations are summa- 



rized in Fig. 3, where we show the results of the scatter- 
ing process at the propagation distance z — 16, similar to 
the corresponding plot in Fig. 2(a). We notice that the 
transverse modulational instability of the stripe starts 
developing for intermediate dark-stripe contrasts, and 
it is accompanied by a subsequent formation of mixed 
(edge-screw) and, later, screw phase dislocations. Con- 
sequently, the stripe breaks-up into vortices of opposite 
charges in a way that resembles the dynamics and forma- 
tion of vortex streets in the hydrodynamics. This phe- 
nomenon can be understood as an effective "unzipping" 
of the dark-soliton stripe by the vortex, that is a mani- 
festation of strongly nonlinear regime in the Aharonov- 
Bohm scattering, a remarkable nonlinear effect. 
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FIG. 3. Snapshots of a dark-soliton stripe scattered by an 
optical vortex soliton, after propagation z = 16 and at dif- 
ferent initial velocities v 2 . The initial offset of the stripe is 
set such that if alone the stripe reaches the origin at z — 8, 
except for the case v = where xo = 2.5. 
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FIG. 4. Vortex pair creation and annihilation during the 
scattering of the stripe (v — 0.1), that gives rise to an effec- 
tively large vortex shift. 

Additionally, we notice that the primary vortex un- 
dergoes a large shift from its initial position after the 
scattering. In order to understand this effect, we have 
studied this process in details (see Fig. 4) and have ob- 
served a novel type of interaction associated with the 
creation/annihilation process of a vortex pair. Indeed, 
during the interaction, a part of the stripe closest to the 
vortex breaks up creating a pair of vortices with oppo- 
site charges [see Fig. 4, z = 4]. Consequently, one of the 
new vortices annihilates with the primary vortex [Fig. 4, 



2 = 5], substituting it by a new vortex of the generated 
pair. This appears as an overall effect of a large vortex 
shift [Fig. 4, z = 16]. 

Conclusions. In this Letter, we have studied an opti- 
cal analog of the Aharonov-Bohm scattering that allows 
illustrating, at a macroscopic level, one of the fundamen- 
tal quantum mechanical phenomena. Moreover, since 
each specific problem on the geometrical phases has its 
own peculiarities, in optics we are able to study, for the 
first time to our knowledge, both weakly and strongly 
nonlinear effects in the vortex-induced wave scattering 
characterized by a significant change in the symmetry of 
the scattered wave and the subsequent stripe instabil- 
ity, unzipping and break-up. Here we have analyzed the 
NLS model in optics since the experimental verification 
of the obtained results is a subject of our following in- 
terest. However the same model is known to appear in 
the macroscopic dynamics of the Bose-Einstein conden- 
sates, so that our analysis applies to that case as well. 
We hope our theoretical results can be further verified 
experimentally, in both nonlinear optics and the physics 
of ultracold condensed gases. 
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